Quantum motion of a neutron in a wave-guide in the gravitational field. 
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We study theoretically the quantum motion of a neutron in a horizontal wave-guide in the gravi- 
, tational field of the Earth. The wave-guide in question is equipped with a mirror below and a rough 

çSJ ■ surface absorber above. We show that such a system acts as a quantum filter, i.e. it effectively 

y—\ ' absorbs quantum states with sufnciently high transversal energy but transmits low-energy states. 

The states transmitted are mainly determined by the potential well formed by the gravitational 
t-H . field of the Earth and the mirror. The formalism developed for quantum motion in an absorbing 

■ wave-guide is applied to the description of the recent experiment on the observation of the quantum 

' states of neutrons in the Earth's gravitational field. 
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I. INTRODUCTION 

Although the solution of the problem of the quantization ofparticle motion in a well formed by a linear potential 
and ideal mirror has been known for a long time Q, H 13- 0, H the experimental observation of such a phenomenon 
K> ' ^ n * ne case °f a gravitational field is an extremely challenging task. 

The elèctric neutrality of neutrons 0, E H El EI is an advantage for this kind of research. Thus, in earlier 
experiments the use of cold neutrons has allowed the gravitationally induced phase-shift of neutrons to be measured 

0111111111111. 

The direct observation of the lowest quantum states of neutrons in the Earth's gravitational field above a mirror 
has recently become possible. The experiment consists of the measurement of the neutron flux through a slit between 
a mirror and an absorber (scatterer) as a function of the slit size. Slit size could be finely adjusted and precisely 
measured. The neutron flux in front of the experimental installation (in Fig^on the left) is uniform over height and 
isotropic over angle. A low-background detector measures the neutron flux at the exit (in FigQ]on the right). The 
main aim of this experiment was to demonstrate, for the first time, the existence of the quantum states of matter in 
a gravitational field. The detailed description of the experiment and a discussion of its reliability and precision can 
be found in refs. [H E d ü EJ |Ü g! H |Ü 

The gravitationally bound quantum states of neutrons and the related experimental techniques provide a unique tool 
for a broad range of investigations in the fundamental physics of partides and fields. These include the Equivalence 
Pri ncip le tests in the quantum domain as well as short-range fundamental forces studies [51 |Ü 0, HI EÜ HI HI 
I30I l3l| and the study of the foundations of quantum mechanics [3^ . l33|. The experiment on neutron gravitational 
quantum states stimulated progress in surface studies (see, for instance. |34l l35j | ) . A short overview of the applications 
can be found in[3(j. 

These studies require clear understanding of the quantum mechanical problem of neutron passage through an 
absorbing wave-guide in the presence of gravitational field. Here we develop a theoretical model of neutron quantum 



FIG. 1: Schematic view of the experiment. From left to the right: the vertical bold lines indicate the upper and lower plates 
of the input collimator (1); the sòlid arrows correspond to classical neutron trajectòries (2) between the input collimator and 
the entry slit between a mirror (3, empty rectangle below) and a scatterer (4, black rectangle above). The dotted horizontal 
arrows illustrate the quantum motion of neutrons above a mirror (5), and the black box represents the neutron detector (6). 



motion in such a wave-guide. 

In Chapter II we summarize the main known facts about a solution of the quantum-mechanical problem for a 
particle in the potential well formed by a linear potential and an ideal horizontal mirror. In Chapter III we discuss 
the main principies of observation of neutron quantum states using the absorbing wave-guide. We show that such a 
wave-guide turns out to be a quantum filter, which absorbs states with high transversal energy and transmits low- 
energy states. These transmitted states are mainly determined by the potential well formed by the gravitational field 
and the mirror. 

The latter condition is a specific feature of our problem, which, to our knowledge, has not been explicitly considered 
in the literature (see, for instance, refs. |33) 13a. l39l I4Ò1 Elj and the references therein, devoted to the theory of the 
interaction of waves with rough surfaces) . Chapter IV is devoted to the passage of neutrons through the wave-guide 
with a flat neutron absorber, as proposed in [3, 13, and Chapter V to their passage with a rough absorber, as proposed 
in |17| . We examine several models for the mechanism of neutron loss as a result of their interaction with an absorber 
and discuss the limits of their validity. 

The final chapter summarizes the conclusions. The results obtained are rather general in character and can be 
applied to different physical problems, involving the transmission of quantum partides through absorbing wave-guides. 



II. QUANTUM BOUNCING ABOVE MIRROR IN THE GRAVITATIONAL FIELD 

Although the results of this section can be found in the handbooks P, 0, 0] it is convenient to have them at our 
disposal here. We start with a well-known problem of a particle bouncing in the gravitational field above a perfect 
reflecting mirror. In the following we consider h — 1. The characteristic for this problem energy scale Eq and length 
scale lo are: 



e = ^?/2 (1) 
lo = ^l/(2m 2 ff ) (2) 

here m represents the particle mass, and g free fall acceleration. In the case of neutrons, which will interest us below, 
these quantities are: 

eo = 0.602 peV, h = 5.871 \im 

The Schròdinger equation, which governs the wave-function of the neutron, confined between mirror and gravita- 
tional field is: 

^72— = X ntp(0 

where dimensionless variable £ is connected with the distance variable z via £ = z/lo, while quantum number A„ 
determines the energy vàlues e n = £oA„. The obvious boundary conditions are: 

ÍPn(0) = (3) 
<Pn(°°) = (4) 

The wave-functions which satisfy the equations above are known to be: 



~ Ai(£- A n ) 



(5) 
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n 




\W KB 


E n , peV 


H n , [irn 


1 


2.338 


2.320 


1.407 


13.726 


2 


4.088 


4.082 


2.461 


24.001 


3 


5.521 


5.517 


3.324 


32.414 


4 


6.787 


6.784 


4.086 


39.846 


5 


7.944 


7.942 


4.782 


46.639 


6 


9.023 


9.021 


5.431 


52.974 


7 


10.040 


10.039 


6.044 


58.945 



TABLE I: Eigenvalues, gravitational energies and classical turning points of neutrons in the earth's gravitational field above a 
mirror 



here Ai is the Airy function (48J. Substitution of J5J) into J3J gives the equation for the eigenvalues A„: 

Ai(-An) - (6) 
The semiclassical (WKB) expression for the eigenvalues is: 

>% KB = {^ (2n-l/2)) ^ (7) 

This approximation gives the eigenvalues with accuracy to a few percent even for the lowest n. 

The asymptotic behavior of the gravitational states' wave-functions in the classically forbidden region £ 3> À n is 
characterized by very fast decay: 

Ai(£ - A„) ~ exp(-2/3(£ - \ n f' 2 ) (8) 

The fast decay of the wave-functions under the gravitational barrier allows us to introduce a well-defined charac- 
teristic distance H n = l X n of a given state, which corresponds to the classical turning point H n — E n /(Mg) of a 
bouncing particle with a given energy. Thus the quantization of energy E n is reflected in spatial distribution of the 
neutron density in the above-mentioned states (hereafter referred to as gravitational states). The scanning of this 
"quantized" spatial distribution of neutron density can be used to observe neutron quantum motion experimentally 
in the gravitational field. 

In TableQJwe present the first seven eigenvalues A„, their WKB approximation X^ KB together with the corre- 
sponding energy vàlues E n = £qA„ and classical turning points H n = lç)\ n - 



III. THE PRINCIPLE FOR OBSERVATION OF THE QUANTUM GRAVITATIONAL STATES 

Here we discuss only the principle of the experimental observation of neutron gravitational states based on the 
concept of neutron tunneling through the gravitational barrier, which separates the classically allowed region and the 
absorber position |20Í |21| . 

A flux of neutrons with horizontal velocity V (from 4 to 10 m/s) was driven through a slit of variable height 
between a perfect horizontal mirror and a highly efficient absorber placed parallel to the mirror. The length L of the 
wave-guide (which varied in different measurements from L — 10 to L = 20 cm) determined the neutron passage time 
T pass _ j^jy ^ 2 10~ 2 s. It was found that when the slit height H was smaller than the height of the first gravitational 
state H\ (see TableUJ the flux of neutrons passing through the slit was indistinguishable from the background. As 
soon as the absorber position was set above H\ a ràpid increase in the flux of neutrons was observed. An analogous 
increase, though less resolved, was observed for the slit heights close to the characteristic state height H2 . This " step- 
like" dependence faded almost completely for higher positions of the absorber, where the flux increased practically 
monotonously. 

We will show here that such behavior of the neutron flux detected at the exit of the wave-guide is what one would 
expect from the qualitative treatment of neutron quantum motion in the gravitational field. In fact, the transversal 
motion of neutrons in the wave-guide can be described as a superposition of the neutron wave-guide transversal modes: 

$(z, t) = J2 C n Mz) exp(-iE n t - T„í/2) 

n 

Here il> n (z) represents the transversal states wave-functions, E n the transversal self-energies and T n the widths of 
these states due to the neutron interaction with an absorber. The neutron flux, detected at the exit of the wave-guide 



4 



is: 

/>oo 

F= \<S>(z,T pass )\ 2 dz 
Jo 

The WKB approach can be proposed for the estimation of the widths of transversal states: 

r„ = p n uj n (9) 

where P n is the probability of absorption of a neutron with energy E n by an absorber during a " one-time collision" , 
while u) n is the frequency of these collisions. The classical expression connecting the frequency of bouncing particle 
and its classical turning point H n is: 




We will use the following simple model for P n . Namely, we will consider P n = 1 when the absorber height H is 
below or equal H n , so that a neutron can "touch" an absorber while it bounces above the mirror in the n-th state. If 
H > H n the probability is equal to the probability of tunneling through the gravitational barrier P — D(E n , H) [Iflj . 
Such a probability has the following form in cases where H 3> H n : 



P n = D(E n , H) ~ exp 



--((H-H n )/l )^ 



(11) 



The spectrum of transversal states depends on the position of absorber H . As long as H > H n the first n states 
can have a long enough lifetime to pass through the wave-guide: 



long 



= J- e xp (i/3(H/l -X n )) 



3/2 



(12) 



The lifetime of all other states with E > E n is approximately equal to the classical time of flight of the particle 
with energy E from the mirror to the absorber. We consider that such a lifetime is short compared to the passage 
time T pass through the wave-guide (which is ensured by the choice of length of the wave-guide L and the horizontal 
flux velocity V) and their contribution to the detected flux is small as far as T short <C T pass . 

The measured neutron flux is: 



N 



F~^|C„| 2 exp(- 



/ ' n 



(13) 



Thus the measured flux exhibits a fast increase when absorber position H is set close to H n due to the exponential 
increase of the n-th state lifetime ljT2*)l . which enables the passage of neutrons in such a state through the wave-guide. 
(A more accurate expression which includes the interference effects between decaying states will be obtained in a later 
section). The expression presented above was used to fit the experimental data: 



JV 

F(H) =Y, A n eM-r pass /r l n ong (H)) (14) 

n=l 

where t^ ™ 9 (íï) are defined by expression (|12fl with X n used as free parameters, while A n = |C„| 2 were used to fit the 
"initial populations" of transversal states. The fitted vàlues of A„ (taking to account the final accuracy of the height 
calibration) are in agreement with the expectation, as given in Table I. The vàlues of A n turned out to be equal, 
except Ai ~ 0-7A n with n > 2. The reason for the approximate equality of the "initial populations" will be discussed 
in the later section. In Fig|21we show the experimental data (2002 year run |20|) and the results of the fit. 

When H ~ Hn and ÍV > 1 a large number of states passes through the wave-guide. This number can be found 
from the WKB expression for eigenvalues Q: 



(15) 



5 




Slit size, (im 



FIG. 2: Neutron count rate at the exit of the wave-guide as a function of absorber position. Circles correspond to the 
experimental data, sòlid line correspond to the theoretical fit. 

Thus for big N 3> 1 the detected flux as a function of H turns to be: 

F(H) ~ {H/l f/ 2 (16) 

The above mentioned WKB expression describes well the flux behavior already for N > 5. The deviation of the 
measured flux from the above expression for small N is due to the quantum character (|14J) of the neutron motion in 
the gravitational field of Earth. Such a deviation (see Fig[2Jl is clearly seen for the first state when quantum formula 
exhibit distinct threshold behavior at H = H\ . However the experimental possibility of resolving higher quantum 
states is restricted by the penetrability of the gravitational barrier. In fact the best resolution of the gravitational 
quantum states is achieved when the flux (|13|l has a step-like dependence on H. This means that the transition 
factor for given state exp(—T pass /T^' n9 (H)) changes from the small value to unity in the range of absorber positions 
H = H n ± S n . The rate of such an increase is limited by the penetration probability through the gravitational barrier 
D(E,H) (nj!. 

For the clear resolution of different quantum states one needs S n <C H n+ i — H„. Under the conditions of our 
experiment 6 ~ 1$ and H n ~ H n . Such an estimation shows that for the highly excited gravitational states (with 
practically N > 5) the difference H n+ i — H n becomes comparable with the uncertainty 8 and thus the step-like 
behavior of the flux is suppressed. To go beyond the above-mentioned qualitative predictions of the resolution of 
gravitational states one needs to take into account details of the interaction of the neutron and the absorber. We will 
return to the discussion of the problem of the resolution of excited gravitational states in the later section. 

IV. FLAT ABSORBER 

The simplest approach in which the properties of the absorber could be taken into account is a model of a flat 
absorber, characterized by the complex Fermi potential. The simplification of the theory in the case of flat absorber 
is due to the fact that, in such cases, motion in a transversal direction is independent of motion in a longitudinal 
direction within the wave-guide. 



6 



A. Passage of the neutron through an absorbing wave-guide 

The Schròdinger equation, which governs the wavc-function <fr(x,z) of the neutron with total energy E passing 
through the wave-guide is: 



i d 2 



i d 2 



2m dx 2 2m dz 2 



mgz + V(H, z) - E 



$(x,z) = 



(17) 



Here x is the longitudinal variable, z is the transversal variable, V(H, z) — Vx(H,z) — Wi(ïl,z) is the complex 
Fermi potential of the absorber dependent on the absorber position H. 

It is convenient to introduce the transversal states ip n (z), which are the eigenstates of the transversal Hamiltonian: 



(e n (H) - iT n {H)/2) 



M*) = 



(18) 



where e n {H) — iT n (H)/2 are the complex energy eigenvalues, dependent on absorber position H. It is worth noting 
that due to the presence of absorption (i.e. the imaginary component of V(z, H))) the above mentioned transversal 
Hamiltonian is no longer self-adjoint. As a consequence, the eigenfunctions ipn(z) are substantially complex and obey 
the bi-orthogonality condition (see |2Í| and references therein) : 



lpk(z)lp n (z)dz = S kn 



(19) 



As one can see the above expression diffcrs from the Standard orthogonality condition in the absence of complex 
conjugation. 

From the qualitative treatment of the previous section one can expect that the lifctime of the neutron in a transversal 
state tp n (z) strongly depends on the absorber position H . The first n lowest states such that H n <C H are weakly 
affected by the absorber and practically coincide with gravitational states ©. Their lifetime is large compared to 
the passage time T pass . The states with H n 3> H are strongly distorted by the absorber. We will show that the 
corresponding lifetimes are short in comparison with r pass and these states totally decay before reaching the detector. 
Consequently, only states with rather small transversal energy and thus small width have a chance of exiting the 
wave-guide. When the absorber position H is reaching one of the characteristic classical turning points H n the 
corresponding state lifctime (and the wave-guide transition factor) undcrgocs fast changcs with H, which allows us 
to monitor this quantum state in the overall flux at the exit of the wave-guide. To calculate the transition factors for 
the given state we expand the two-dimensional wave-function <fr(x, z) in the set of basis functions ip n (z): 



(20) 



The functions Xn (x) play the role of longitudinal wave- functions of neutrons in transversal state n and can be found 
by substitution of (|20|l into the Schròdinger equation <|17D with the use of (|19fl : 



1 d 2 
2m dx 2 



e n (H)-iT n (H)/2-E 



Xn(x) = 



(21) 



The solutions of l|21() corresponding to the quasi-free longitudinal motion of neutrons in the wave-guide are: 

Xn{x) ~ exp(ip„x) (22) 



Here p n — ^/2m(E — e n (H) + íT n (H)/2) is the complex longitudinal momentum. 

As we have already mentioned, only states with small transversal energy can reach the detector. In our case the 
full energy is much greater than the transversal energies: 



|e n (H)+<r n (ff)/2)| «£ 



Thus we can write for momentum p n : 



~ e n (H)-iT n (H)/2 ^ e n (H) ~ iT n (H)/2 
Pn ^ V2mE(í — ) = P 
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where P = \/2mE and V = P/m. Due to the positive imaginary part of p n , each of the longitudinal wave-function 
decays exponentially with x inside the wave-guide: 

, v - g „(g) 3 ; + »r„(g)x/2 ^ 

X«W ~ exp(z — ) exp{iPx) 

Taking into account that r pass = L/V, we obtain for the wave-function 9(x = L, z) at the exit of the wave-guide: 

*(as = L,z) = cxp( ï PL)^C„exp(-ze Il (ií)TP ass )exp(-í^ffi^)^(z) (23) 

n 

Here C n are expansion factors, determined by the particular form of the wave-function at the entrance of the 
wave-guide: 

poc 

C n = V(x = 0,z)i; n {z)dz (24) 
Jo 

B. Expansion factors 

It is worth mentioning that strictly speaking |C„| 2 cannot be interpreted as the initial population of the certain state 
ipn(z)- In fact, as long as the Standard orthogonality condition is not vàlid for the eigenfunctions of the not-self-adjoint 
Hamiltonian (ip n \ipk) ^nk we have: 

/>oo 

V|C„| 2 ^ / \9(x = 0,z)\ 2 dz 

Consequently, to find the measured neutron flux at the exit of the wave-guide one has to take into account that: 
F = / |*(a: = L, z)\ 2 = ]T C^C^n^u) eM~^k ~ e„)r pass ) exp(- 1 " + k) ) ± ^ \C n \ 2 exp(~r„r^ s ) 

, ' n,k n 

(25) 

The appearance of the interference terms C*C/ c ( - í/vi|V'fc) is not surprising. In fact, the states \ijjk) ^ve not stationary 
states with certain energy. Due to final decay width these states are in fact time-dependent and can be expressed as 
superpositions of stationary states with certain energy. The contribution of the mentioned interference terms to the 
flux can be interpreted as oscillating in time transitions with frequency u)^k — e n ~ £ k between the true stationary 
states. 

However for the observation of the interference terms above, a rather narrow distribution of neutrons is required in 
the longitudinal velocity. Should such longitudinal velocity distribution be broad, the interference terms are canceled 
after averaging over such a distribution and the " Standard" expression for the flux is restored: 

F = J2\Cn\ 2 eM-^nr pass ) (26) 

n 

Indeed, the interference terms C*Cfc appear in the expression for the measured flux 1)25(1 multiplied by 
exp(zcj„fcr ï ' a;is ). If the initial flux has distribution f(V) in the longitudinal velocity, the contribution of the interference 
terms averaged over such a distribution would be: 

C* n C k (ip n \ip k )eM^nkT pass )eM- — * )f(V)dV 



In the case broad velocity distributions, such that: 

AV , 

\ r^pass \ — 1 
-— > (T F CJnk) 

the contribution of the interference terms is canceled due to the fast oscillating term exp(io; íl fcT pass ). To observe 
the interference contribution between the first and second states, the velocity resolution in the conditions of our 
experiment should be better than 10%. This limitation is less severe for excited states. 
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Let us now turn to the problem of the initial " population" of the gravitational states where the initial flux has 
broad distribution in transversal momentum. In such a case [30l | the modulus square of expansion cocflicient |C„| 2 
can be found from the following equation: 

|C„| 2 = í(ip n \k)(k\^ n )exp(-k 2 /k 2 )dk (27) 

where ko is a characteristic width of the transversal momentum distribution and we have used "bra-ket" notation for 
the matrix element (VvJfc) == f ip n (x) exy(ikx)dx. 

It has been shown in |3fl|3Ï],|5(j that if k l ^> 1 the squares of the amplitudes of the lowest states are practically 
equal: 

|Cn| 2 ~ í-o(rV) 

k t 

In the conditions of our experiment the corresponding value kalo ~ 50 and thus the approximation of a unificd 
population of lowest states is well justified. 

Indeed, having in mind fast oscillations of the integrand in l|27|) |C„| 2 becomes very small if k > k c , where k c ~ 1/H n 
is the characteristic momentum of the gravitational state with spatial extension H n . As far as the distribution over 
k in the initial flux is practically uniform for k < l/H n <§C ko the expression (|27f) can be rewritten as: 

\C n \ 2 = J \^ n \k){k\4>n)dk = (V„|Vn) = 1 

and we return to the statement of the uniform distribution. It is worth mentioning that the same averaging over the 
initial transversal momentum distribution applied to an evaluation of the interference term C*Ck gives: 

c*c k = (vg^) 

As we have already mentioned, this matrix element is nonzero for those states which are affected by the absorber 
and depends on absorber position H. Given the above arguments we can rewrite the expression for the measured flux 
as a function of H (I25() after averaging over the transversal momentum of the initial flux as: 

F(H) = J2cM-rn(H)T» ass )+ Yl 2 Re (^ n \A) 2 exp^^g)^)) exp(- (F " (g) + ^ígljlüü ) ( 28 ) 

n—l n,k>n 

In the case of a broad longitudinal velocity distribution in the incoming flux, only the first term in this expression 
is important. 



C. Transition factor 



Once the expression (I28f) has been obtained, the problem of calculating the neutron flux at the detector position is 
transformed into the problem of calculating the eigen-energies e„ and their widths T n (H) of transversal states as a 
function of absorber position H. 

The realistic Fermi potential V(H, z) of the absorber material is characterized by the depth of order 10 _8 eV i.e 
much greater than the characteristic energy ÍO -12 ^ of the lowest gravitational states. The diffusion radius p of such 
a potential, i.e. the distance where the strength of potential rises from zero value in the free space to its final value 
inside the media, is much less than the characteristic gravitational wave-length Iq. In such a case the properties of 
the absorber can be precisely described by one parameter, namely the complex scattering length a, whose imaginary 
part accounts for the loss of neutrons due to absorption. (We use hereafter the following definition of the scattering 
length a — Umfe_»o(l — S)/ (2i/c), where k is neutron momentum and S is the reflected wave amplitude). An analytical 
equation for the eigen-energies of neutrons bouncing in the gravitational field between mirror and absorber which is 
positioned at distance H above the mirror can be derived. We refer the reader to Appendix A for the details and 
present here the final expression for the eigenvalues A„ (H) : 

Ai(-À n ) = Ai{H/lo - A n ) - à/l Ai' (H/lp ~ K) . . 

Bi(-A») Bi{H/l ü - A„) -à/l Bi'(H/l - A„) 1 ' 

Here à = a — H plays the role of "the scattering length on the diffuse tail" of the potential. Let us note that this 
expression is vàlid for any Fermi potential with a small diffuse radius, as long as the corresponding scattering length 
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a is small compared to the gravitational wave-length lo. If the position of absorber H ^> H n the right-hand side of 
the equation (|29[) becomes exponentially small: 

'0/ 

One can easily recognize in the right-hand side exponent the penetration probability through the gravitational 
barrier. Taking into account the smallness of such a probability one can find the correction to the gravitational 
eigenvalue A„ due to the small, but nonzero, possibility of penetration under the gravitational barrier to the absorber. 



Ai(-A n ) _ 1 
Bi(-A„) 



2 6XP 



--((H-H n )/l )^ 



Bi(-A°) 

AA n = / ; „ , exp 

2Ai'(-A° 



--((H~H n )/l Q ) 3 / 2 



2 v / (H-H n )/l Q - + l 



(30) 



where A^ are unperturbed eigenvalues determined by ijïïj). The width of the n-th state due to penetration under the 
gravitational barrier and absorption turns out to be: 



T„ ~ 2 l -^eoJ^-V(H-H n )/l e X p 

'0 V tln 



- ï {{H-H n )/lo) 3/ * 



(31) 



We have used in the derivation of this expression the semiclassical approximation for the Airy function. The physical 
sense of this expression for the decay rate becomes clear after comparison with the semiclassical expressions (J§J. 
Taking into account the expression Í|1C)|) for the classical frequency uj n we can rewrite ()31(l as: 



4Í^ Wr V(iï-iïn)A)exp 
to 



; ({H-H n )/l f' 2 



The neutrons penetrate through the gravitational barrier into the absorber, the corresponding probability (|ll|l is 
exponentially small. This probability is multiplied by the classical bouncing frequency in given state n. The properties 
of the absorber itself appear in the above expression through the ratio 4 I " la ^ yj (H — H n )/lo. Later we will show that 
it coincides with a general expression for the absorption probability of slow quantum partides on the short-range 
absorbing potential. Thus the intuitive formula 10 we used before is justified. 

One can introduce the characteristic absorption time: 



lo 



'Hn 



2| Imajeo V 



For efficient absorption one needs T P ass / T ^ bs ^> i j which puts the following requirement for the imaginary part of 
the scattering length Ima: 



Imal > 



lo 



2e rP as 

In the conditions of our experiment (r pass = 2 x 10 -2 s) the above requirement means that 

Ima ^> 0.05Zo ~ O.Sfim 



(32) 



(Let us mention that at the same time the scattering length approximation used above is vàlid only for vàlues of 
|5| < lo) 

This treatment shows that for the clear resolution of quantum states the most favorable absorbers are those with 
the largest possible scattering length. 

Let us see how the scattering length discussed above is connected to the properties of the absorber's Fermi potential, 
namely its complex depth and the diffusion radius p. We will study the case of the complex potential of the Woods- 
Saxon type: 



V(z,H) 



U exp(— icp) 



1 + exp((H - z)/p) 

In the following U > 0. It can be shown that the scattering length on such a potential is given by: 

r(l + pn) 



a = H - - + 2p[j 



(33) 



(34) 



10 




0^ 1 r~'ï~"T ï , i , 1 , 1 , 1 , 1 

10 15 20 25 30 35 40 45 

H, (nm) 



FIG. 3: Lifetime of the first three gravitational states as a function of absorber position in a typical arrangement of our 
experiment 



Here 7 0.577 is the Euler constant, k = exp(— i(p/2)y/2mU and T(x) is the Gamma-function. 

An important limiting case is the case of the deep complex Fermi-potential, namely p\n\ 3> 1 and p\ Iuík\ 3> 1. It 
follows from l|34l) that in such a case: 

Ima = — ptp (35) 

and is independent of the depth U of the complex Fermi potential. It has been shown in |5l| that such behavior of the 
imaginary part of the scattering length is universal for deep complex potentials with the exponential tail. For such a 
strong absorbing Fermi potential the neutron is completely absorbed on the tail of the complex Fermi potential; the 
properties of the inner part of the absorber therefore lose their importance. The only way to increase the scattering 
length in such a limit is to increase the diffuscncss p. 

In an other limiting case, when the diffuseness is so small that pn <C 1 the scattering length becomes: 

1 7T 2 , 
Im a — — Im —p k 

k ò 

The leading term in the above expression is —1/k. One can check that it coincides with the scattering length, 
characterizing the low-energy reflection from the step-like potential t/exp(— iíp)Q(z). Indeed, in the case of pn <C 1 
the neutron can penetrate through the narrow exponential tail of the complex potential into its core without significant 
losses. One can see from the above expression that for weak absorbers with depth U ~ £0 the imaginary part of the 
scattering length becomes as large as Iq. 

It is important to mention here that the absorption of the ultra-cold neutrons by the complex potential is closely 
related to the so-called quantum reflection 52, 53] of ultra-slow neutrons from the fast changing complex Fermi 
potential. The reflection probability R in case of slow quantum partides |H^ | impinging on the absorber at normal 
incidence with momentum k can be written as follows: 

fí=l-4fc|Ima| (36) 

while the absorption probability P is: 

P=l- R = Ak\lma\ (37) 

The smaller the p (and | Ima|) the better the reflection and the weaker the absorption. One can see that the limit 
of up — > 0, p/k — > corresponds to the case when the absorber is replaced by an absolutely reflecting mirror. The 
above mentioned high reflectivity of a fast changing potential is a general quantum mechanical property. 

The numerical calculations verify the above conclusions. 

The vàlues U, p and <j) of potential (|33|l were chosen to be U = 10~ 8 cV, p = 1 pm and <f> = 3n/4 (this corresponds 
to attractive potential with absorption). 

In FigEl we plot the lifetimes r„ of the first 3 states as a function of absorber position H . In Fig0| we show the 
corresponding evolution of the real part of energy e n (H). One can see the fast increase in the lifetimes at certain 
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FIG. 4: Energy of the first three gravitational states as a function of absorber position in an arrangement typical for our 
experiment 
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FIG. 5: The relative neutron flux as a function of absorber position for different absorber difuseness in an arrangement typical 
for our experiment. F* is the flux calculated at absorber position H — 45 pm and diffuseness p — 1 p,m. 



vàlues of H, close to H ni in agreement with the qualitative predictions of the previous section. The real part of the 
energy quickly approaches its limiting value equal to the energy of the gravitational state when H > H n . 

As long as these states are treated separately the clear evidence of the fast changes in the lifetime, as a function of 
absorber position can be seen. However, the overall plot of flux intensity FigEl where all these states are taken into 
account simultaneously shows that the step-like dependence is suppressed, except for the first step at H = H\ and 
partially for the second step. 

To achieve much higher absorber efficiency the diffuse radius p should be significantly increased. Another way is to 
reduce the depth of absorber Fermi potential to the level of 10 -12 eV, the characteristic scale of gravitational states 
energies. 

Absorbers with optimal parameters can be obtained if their surface is corrugated. In fact such an absorber was 
used in the experimental set-up. The zone of such a corrugation can be considered as a low density media with an 
extended diffuse radius of Fermi potential. In the following we will study the neutron passage through an absorber 
with a rough surface. We will show, however, that the main loss mechanism in such a case is due to non-specular 
reflections from the rough edges of absorber. 
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D. Zero gravity experiment 

We will study here the important case of the neutron passage through the wave-guide formed by the mirror and 
the absorber in the absence of the gravitational field. The case is interesting from two points of view. On the one 
hand, a comparison of the transition factors with and without gravity clearly shows the role of the latter EU • O n 
the other hand the "zero" gravity experiment (which simply means the installation of the mirror and the absorber 
paral·lel to the gravitational field) enables independent measurement of the mirror and absorber properties. 

Let us first mention that the neutrons' motion transversal to the direction of the mirror (and the absorber) is 
quantizcd. Neutron states of this type, localized between the mirror and the absorber, will be referred to as "box- 
like". However due to the loss of neutrons inside the absorber such states are no longer stationary states; they are 
quasi-bound states with finite life-times (width) . The existence of quasi-bound states in the presence of an absorber 
is a consequence of the phenomenon mentioned above as quantum reflection from the fast changing absorber Fermi- 
potential. In fact the partial reflection of neutron waves from the absorbing potential leads to the formation of the 
standing wave (i.e. quasi-bound state). The more efficient the absorber the smaller the amplitude of the reflected 
wave and the shorter the neutron life-time. In the case of full absorption of the neutron wave (which means that the 
amplitude of the reflected wave is exactly zero) no quasi-bound state can exist. 

With these remarks we can now turn to the calculation of the neutron flux through the wave-guide: 

F = ^|C n | 2 exp(-r n r í,aS5 ) (38) 

n 

Here n is the quantum number of the quasi-bound box-like state. In the above expression we neglect for the 
moment the contribution of the interference terms. As we have shown, this is possible when the longitudinal velocity 
distribution is rather wide. In the following we will also assume a wide distribution of the incident flux over transversal 
momentum (orthogonal to the mirror and absorber). We have already established that in such a case the first 
states are populated homogenously. The number Nh of homogenously populated states can be estimated from the 
condition that the characteristic momentum of the box-like state k c ~ n/H is equal to the spread of the transversal 
momentum distribution fco in the incident flux: 

N h ~ k H (39) 

Let us now turn to the calculation of the widths of certain neutron states, confined between mirror and absorber. 
As in the case of the gravitation states we assume that the absorber Fermi potential can be characterized by a complex 
scattering length a, which is possible when k n p <C 1, (where k n — \/2mE n is the neutron momentum in given box-like 
state with energy E n ). To obtain the complex energies of the box-like states we note that the neutron wave function 
in the region where the absorber Fermi potential can be neglected is: 

* b (z) ~ ún{k n z) 

Such a wave function can be matched with the asymptotic form of the neutron wave-function inside the absorber 
at distances H — l/k n <C z <C H — p, where absorber potential vanishes. The general asymptotic form of the 
wave-function in this region is: 

H-z 

f„(z)~l + ^ 
a 

where 5 = a — H is the "diffuse tail" scattering length. The matching of the wave-function and its derivative leads to: 

K = (40) 
H — a 

ir 2 n 2 2-K 2 n 2 Reà 
n ~ 2nHP + 2miJ3 ( } 

„ „ llmal 7r 2 n 2 |Ima| 

r„ « AEj l -=A l - (42) 

We will show that the dependence l|42|l will play a crucial role in establishing the wave-guide transition factor 
dependence on H. Let us note here that the expression for the width of the box-like state (|42H is a consequence of the 
quantum reflection from the fast changing tail of the absorber Fermi potential. To see how the quantum reflection 
phenomenon is connected to the width of the box-like state let us return to the semiclassical expression for the loss 
rate: 

r~ Lü n P 
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where w„ is the clàssic frequency of collisions with the absorbing wall and P is the probability of absorption in a " one 
touch" collision. The expression for the collision frequency with one of two walls is 



2H 2mH 2mH 2 
The probability of absorption P (|3"T|) turns to be: 

„ „ , , 47rn|Ima| 

P = 1- R = 4k n \Ima\ = l - l - (43) 

Combining the above results for the frequency ui n and absorption probability P we return to the expression 1(42(1 . The 
quantum properties of neutron motion appear here through the energy dependence of the absorption probability 1431) 
and quantization of the box-state energy (momentum) . 

Integrating the results for C„ and T n into the expression for the flux ((26(1 we obtain: 

n v ' 

Here Fq is the normalization constant, characterizing the intensity of initial flux. 

One can see that the number of states passed through the wave-guide is obtained from the condition: 

r„r pass ~ i 

which gives 



27r v /|Ima|Tï> ass 

Hereafter we expect that the number of homogeneously populated states Nh 139(1 is greater than N pass : 

k H > N pass 

From the expression ((44(1 it follows that the wave-guide transition coefncient is determined by the characteristic 
absorption constant: 

7r 2 r pa SS | Ima | 
£ = 4- 



2mH 3 

which is connected with the number of states passed through the wave-guide via: 

N pass 

There are two important limiting cases. 

The first case, which we call the "strong absorption" limit £ > 1, means that a maximum of one state only can pass 
through the wave-guide, i.e.: 

N pass _ 1/^/1 < l 

In this case the neutron flux is a rapidly increasing function of H: 

The opposite case, which we call the "weak absorption" limit £ 1 means that a large number of states can pass 
through the absorber: 

N pass = l/^ > 1 

In this case the summation in expression ((4411 can be substituted by integration, which gives: 

F^Fq (46) 
2^/2^1 lma\TP ass 



14 




FIG. 6: The relative neutron flux in the presence of, and in the absence of, gravity. F* is the flux with gravity calculated for 
absorber position H = 45 fim, | Ima| = 2 /im and T paas = 0.02 s. 



It is worth mentioning that JT 3 / 2 -dependence is a consequence of the quantum threshold behavior that determines 
the energy dependence of the absorption probability of ultra-cold neutrons (|37() . This expression is vàlid in the so- 
called anti-classical limit k n à <C 1. In the opposite case, when k n à > 1 the absorption probability energy dependence 
differs from (|43|l . In particular, if the absorption occurs with unit probability for each collision P ~ 1 it is easy to 
establish: 



T„ w nn/(2H) 

The substitution of this expression into (|26J) results in the H 2 dependence of the flux instead of iï 3 / 2 . Note that the 
large value of the absorption probability of ultra-slow neutrons can be achieved only if there is a large imaginary part 
of the scattering length k n \ Ima| ^> 1. 

Let us also note that we restrict ourselves with the condition of a homogeneous population of box-like states 
kç)H > N pass , so far £ cannot be smaller than: 

Uin = ^2 

Obviously, in the limit of very small £ <C £min , when absorption can be fully neglected, the flux passed through the 
slit starts to be proportional to the slit size: 



F~ H 



which means that all the neutrons that enter the slit pass through it without losses. So far, depending on the efficiency 
of the absorber one can get different flux dependence on the slit size H . 

A comparison with the gravitational case results in the following conclusions. 

First, for small slit sizes, both flux curves, seen as a function of H manifest fast increases in the vicinity of the 
characteristic value H c . However in the case of gravitational states this critical slit size is determined by the "height" 
of the ground gravitational state H c ~ H\ , while in the case of zero gravity it is fully determined by the properties of 
the absorber, namely the imaginary part of the scattering length | Ima| and the passage time T pass : 

H c = (27r 2 r pass |Ima|/m) 1/3 

Secondly, in the presence of gravitation the flux exhibits a step-like dependence on H with increasing slit size, and 
tends to H z / 2 dependence for large H. Such step-like behavior is more pronounced for larger diffuseness of the absorber 
(larger | Ima|) or for longer passage time. In case of zero gravity the flux increases with H monotonously. Such an 
increase has power law dependence in the limit of large H . Depending on the absorber efficiency the corresponding 
exponent can vary from 2 (full absorption) to 1 (no absorption). 

We plot on Fi^üjthe neutron flux in the presence of and in the absence of gravity for | Im a\ =2 fim and T pass = 0.02 

s. 
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E. Inverse geometry experiment 

Another way to clarify the gravitational effects and to measure the efficiency of the absorber is to exchange the 
position of the mirror and the absorber in the experimental setup. Here we will study such an inverse geometry 
experiment, in which the absorber is placed below and the mirror above. 

First we will study the modification of the gravitational energy vàlues due to the interaction with an "absorbing 
mirror" . As we have already shown, as long as the distance p where absorption takes place is much smaller than the 
gravitational wave-length h, such an interaction can be characterized by only one parameter, namely the scattering 
length a -C h, regardless of certain details of the absorber Fermi potential. 

The modification of the eigenvalues A n due to the interaction with an "absorbing mirror" can be obtained by 
matching the wave-function of the neutron, reflected from the absorber, which large z asymptotic form (z 3> a) in 
the case of small neutron energies can be written as follows: 

ip(z) ~ 1 — z/a 

with the gravitational wave-function Ai(z/l — A), where A is a modified eigenvalue. We take into account that in the 
matching region z/lo 1 we obtain the following equation for A: 

^ = -a/h (47) 
Ai'(-A) 

As far as |a/Zo| *C 1 we get the following expression for the modified eigenvalues A„ accurate up to the first order 
of small parameter |a/Zo| : 

À„ = A„ + a/l (48) 

From the above equation we obtain the following modified energy levels: 

E n = e(X n +Rea/h) (49) 

T n = 2e— — ^ = 2m 5 |Ima| (50) 
h 



If we use the expression (|35|) for the scattering length on the deep [2p\j2mU ^> 1) imaginary (ip — tt/2) exponential 
potential, we obtain for the width of the gravitational state: 

T mv = mgnp (51) 

One should mention that the width of the gravitational state l|50|l is independent of the energy (for such states 
that \J2mE n p <C 1). This can be easily explained by the following simple arguments. The frequency of the neutron 
bouncing above the surface in the gravitational field is lj ~ 1/yE, while the probability of the absorption P = 
4fc|Ima| ~ VE. Combining these two variables we get the energy-independent expression for the width T = luP. 
This means that all the gravitational states which are not affected by the upper mirror (H n <C H) decay at the same 
rate l|51fl . The corresponding lifetime in case of p = 1 pm is 



t — 1.7 x 1CT 3 s 

which is much smaller than the passage time T pass = 0.02 s. The expression H50(l manifests the very important 
property of the neutron bouncing above the " absorbing mirror" , namely the factorization of gravitational properties, 
which appears through factor mg and the absorber properties, characterized by | Ima|. 

Let us now understand the behavior of the transversal states with much higher energy E ^s> mgH . For such high 
energies the influence of the gravitational field can be neglected (for neutron motion between mirror and absorber). 
The corresponding states can be treated as the previously studied "box-like" states of the free neutron, confmcd 
between the absorber and the mirror; their widths are given by l|42() . 

As long as we study the transversal states with E n mgH their lifetimes are much smaller than those of the 
gravitational states already considered, and their contribution to the neutron flux at the exit of the wave-guide can 
be neglected. 
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FIG. 7: The relative neutron flux in the direct and inverse geometry experiment. F* is the flux in the direct geometry case, 
calculated for absorber position H = 45 /jm, | Ima| = 2 /im and T pass = 0.02 s. 

The two limiting cases studied above naturally follow from the equation for the eigenvalues for the inverse geometry 
experiment (see Appendix A for details of the derivation): 



a [M(H/l - A„) Bi'(-A„) - Ai'(-A„) Bi(H/l - A„)] - Ai(-A„) Bi{H/l - A„) - Bi(-A„) Ai(H/l - A„) (52) 



We come to the conclusion, that, with mirror position H H%, the measured neutron flux is mainly determined 
by the gravitational states passed through the wave-guide such that H n < H . The number of such states is given by 
(|15fl and thus the dependence of the flux on H is given by l|16|) . The ratio of the fluxos in the "direct" and in the 
inverse geometry experiment turns to be: 



The results of comparison of neutron fluxes in direct and inverse geometry experiment is shown in Fig0 

This difference in the fluxes clearly shows the role of gravitation in the passage of the neutrons through the wave- 
guide. On the other hand it enables us to measure the efficiency of the absorber. It is also interesting to note 
that if |Ima| is known by independent measurement (e.g. from the zero gravity experiment discussed above), the 
measurement of the lifetime of neutrons bouncing in the low gravitational states above an absorbing surface will give 
direct access to the gravitational mass of neutron m and will allow us to apply the quantum equivalence principle 
test. 

The inverse geometry measurements were performed during one of the first runs of neutron gravitational states 
experiment |l9j. The obtained results verify strong suppression of the flux in inverse geometry experiment case in 
agreement with (|53l . 



The previous analysis shows that in order to increase the efficiency of flat absorbers one needs to use substances 
either with a Fermi potential of large diffuse radius or of very small depth (U ~ 10~ 12 eV). The construction of such 
absorbing materials is rather problemàtic. An alternative way to increase absorber efficiency is to use an absorber 
with a rough surface. In the wave-guide experiments an absorber with a rough surface was used with a roughness 
amplitude of about 2 [im. In this section we will study the role of roughness in the neutron loss mechanism. 



F mv /F dlr ~ exp(-2m 5 |Imo|r pas5 ) 



(53) 



V. 



ROUGH SURFACE ABSORBERS 



A. 



Effective potential approach 



The rigorous study of neutron interactions with a rough surface requires solving the two-dimcnsional problem, where 
the neutron-surface interaction is described by a rather complicated function V(x, z). The radical simplification of such 
a problem is possible via introduction of effective one-dimensional potential V e ff (z) 50]. The simplest assumption 



17 



enabling us to calculate such a potential is the following. We expect that the longitudinal kinetic energy of neutrons 
Pq/2M is sufficiently superior to the characteristic value of the Fermi potential V(x, z) of rough edges. The first order 
Born correction to the longitudinal kinetic energy of neutrons due to the interaction with the rough edges would then 
be: 



AE(z) = -J V{x,z)dx (54) 

where the " normalization length" L is selected to be much greater than characteristic correlation length of roughness. 
This correction to the longitudinal energy, being a function of z, plays the role of effective potential V e ff(z) = AE(z) 
in the equation for the neutron transversal motion: 

1 d 2 

-2^^{z) + AE(z)<p(z) = {E- P 2 /2m)<p(z) 

The physical meaning of expression (|54|l is transparent; it is the potential of media with reduced density. In 
particular if one models the roughness by the periòdic gratings with z dependent width d(z) and period L, than the 
effective Fermi potential is V e ff(z) = Ud(z)/L, where U is the corresponding Fermi potential of flat surface. The 
benefit of this approach is the ability to connect the one-dimensional effective Fermi potential with averaged shape 
properties of roughness and realistic Fermi potential of absorber substance. We will not take this case any further, 
since the main results have already been discussed in the section devoted to the flat absorber. 

The above approximative model can be justified for the longitudinal energies of neutrons much higher than the 
Fermi potential of rough edges. A very important effect, which is not taken into account in this simplified approach is 
the possibility of non-specular reflections, i.e. the energy exchange between the horizontal and vertical motion of the 
neutrons. (They appear in the second order Born approximation) . In the following we develop the non-perturbative 
formalism in which such effects would be taken into account. 



B. A time-dependent model for the neutron loss mechanism 



In the previous analysis we found that only those neutrons which have sufficiently small transversal energy do not 
penetrate through the gravitational barrier into the absorber and thus are not absorbed in the wave-guide. The role 
of the absorber's roughness is to transfer a significant portion of longitudinal energy into transversal energy during 
non-specular reflection from the rough edges. Thus the neutron interaction with the rough surface absorber results in 
mixing of states with diffcrcnt transversal energies. As long as the states with large transversal energy have very small 
lifctimes, such a mixing results in a loss of neutrons. Here we will study this loss mechanism within the time-dependent 
model. 

We will study the neutron passage through the wave-guide in the frame, moving with horizontal velocity V of the 
incoming flux (we suppose that this velocity is well defined) . The rough edges of the absorber surface can then be 
treated as a time-dependent variation of the flat absorber position. This means that the neutron loss mechanism in 
such a model is equivalent to the ionization of a particle, initially confined in a well with an oscillating wall. 

The time-dependent Schròdinger equation for the neutron wave-function is: 



z) 



2m dz 



^^IÏ + m 9Z + V{z,H{i)) 



$(t,z) 



The time-dependence appears here through the time-dependence of the absorber position H(t). 
The boundary conditions are: 

${t,z = 0) = 
$(í,z = oo) = 

It would be convenient here to introduce time-dependent basis functions: 



(55) 



(56) 
(57) 



(j>n(t,z) = ip n (H(t),z)exp(-- / e n (H(T)dr) 



where ip n (H,z) and e n (H) are complex eigenfunctions and eigenvalues of transversal Hamiltonian (|18(l with fixed 
absorber position H. 
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The total wave-function <ï>(i, z) can be expanded in the set of functions: 

$(t,z) = y2c n (t)ÍJ n (H(t),z)eM-i f e n (H(r))dr) (58) 

The equation system for the expansion factor C n (t) is: 

dC n (t*) dH v "\ r , . . 

— t: — = — — ) C k a nk exp[-iu; nk {t)\ (59) 

«nfc = -«fen = y Vn{H,z) ^ (60) 

where 



Wnfc(*)= / [e fc (-ff(r)) - £ n (íí(r))]rfT 
Jo 

Note that the derivation of equations (|59|l and (|60() requires the biorthogonal condition (|19fl . 

The initial conditions Cfe(0) are determined by the overlapping of the incoming flux with the basis functions: 



C7 fc (0)= / ${t = 0,z)i> k (H,z)dz 
Jo 

The solution to equation l|59(l together with the above initial conditions enables us to obtain the wave-function <ï>(í, z) 
at t = r pass and to calculate the measured flux: 



F 



/ \<S>(T PaS8 ,z)\ 2 dz 

Jo 



The equation system (|59|l can be very much simplified under the following assumptions. 
First, let us suggest that the absorber position time-dependence is harmònic: 

H(t) = H + fesin(wí) 

where b is the roughness amplitude and frequency u> = V/d, with d being the spatial period of roughness. It is 
known that in such cases it is only the states in the equation system l|59|l obeying the "resonance" condition: 

\e k -s n \~uj (61) 

which are effectively coupled. As long as the transversal states have the widths this resonance can not be exact. 
However we will restrict our treatment to only two coupled states. 

Second, we expect the roughness amplitude to be so small that the following approximation is vàlid: 

~ fecosM dH \ H=Ho 

Thirdly, we will consider that to S> \s n \, so that the low lying gravitational state ij) n is coupled with the very highly 
excited state with energy Ree^ ^> MgH. The gravitational potential can be neglected in comparison with such high 
energy; we are thus dealing with a "box-like" state. Its energy and width is given by (|41|l and (|42J) . As the width 
of such excited states is much bigger than the width of the low lying gravitational state ip n we can neglect the latter 
and suppose that neutrons in the low-lying gravitational state are elastically reflected both from the mirror and the 
absorber. This results in the following boundary conditions for the gravitation state wave-function ip n (H, z): 

^ n (H,z = 0)=ip n (H,z = H) = 

The eigenfunction ip n (H,z) and the eigenvalue e n (H) — So\ n (H), determined by the above boundary condition 
are: 

MH,z) ~ Bi{-X n {H))M{z-\ n {H))-M{-X n {H))Bi{z-X n {H)) (62) 
Ai(H/lo-X n (H))Bi(-X n (H)) = Ai(-X n (H))Bi(H/l -X n (H)) (63) 
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Finally we come to the equation system with only two coupled equations: 



C (t) = -ibwd^a^expKw-woOí] (m 
C?i(t) = |6wC (í)a(ir)exp[-<(w-woi)t] 1 ' 

with 

= Jo dïï — 

and u>oi = E\ — Eq. 

In the above expressions index labels the low lying gravitational state, while index 1 labels the excited fast 
decaying box-like state with complex energy E\ = Rei?i — iT /2. 

A very convenient expression |55j can be obtained for coupling matrix element a(H) (see Appendix B), namely: 



A — Ai 

The benefit of such a simplified equation system is that it enables an analytical solution. Taking into account initial 
conditions Co(0) = 1 and Ci(0) =0we get: 

C (í) = exp(-Ç) (cos( 7 í/2) + ^sin( 7 t/2)) (66) 

C x (í) = -i^í^expfÇjsm( 7 í/2) (67) 



Here T is the width of the "box-like" state and 7 = l/2y/b 2 ui 2 a 2 (H) — T 2 . (Note the exponential increase of Ci(t). 
This does not yield in nonphysical result, as far as in the expression for the wave-function (|58jl C\{t) is multiplied 
by decaying exponent exp(— iE\t). However, as we mentioned before, | C (í) | 2 and |Ci(í)| 2 cannot be interpreted as 
probabilities to find a system in certain quantum state). 

We are interested in the evolution of the gravitational state. Two important limiting cases are: 

|C (Í)| 2 - exp(-^), if n 2 /T 2 « 1 (68) 



|C* (í)| 2 -> cxp ( ) cos 2 {üt/2 ~ip)/ cos(^), ïïn' 2 /T 2 > 1 (69) 



4r 
n 

~2 

Here tp = arctan (r/(2 7 )) and Q 2 = b 2 uj 2 a 2 (H) 

The quantity Çl plays the role of "transition frequency" between two states. It is proportional to the roughness 
amplitude b and depends on the averaged absorber position H via the coupling a(H). The coupling a(H) decays 
rapidly as soon as H > H n , where H n is the classical turning point for the low-lying gravitational state. 

When r ^> fi the decay rate r n of the n-th gravitational state , according to 168(1 . is r„ = íí 2 /(4r). Using the 
asymptotic expressions for a(H) (see Appendix B) one can get the following expression for the decay rate in case 
H » H n : 



lo b 2 H — H n 



i((7í-ií„)//o) 3/2 



(70) 



where Ima is the imaginary part of the scattering length of the neutrons on the flat absorber Fermi potential. 
The expression (|7UH should be compared with the analogous formula for flat absorbers 131fl . One can see that 



«off 



16 Ima 



plays the role of the effective scattering length of the rough surface absorber, which is proportional to the square of 
the roughness amplitude. 
The time 

abs 8/ |Ima| 
b 2 s 
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FIG. 8: The relative neutron flux as a function of the slit height in the time-dependent model. F* is the flux calculated for 
absorber position H = 40 /im, 6=1 /im, | Ima| = 0.1 fim and T pass = 0.02 s. 




plays the role of the characteristic absorption time in our problem. 

It should be noted that the above results are true for H > H n and "weak" coupling. When the absorber position 
H < H n the coupling is large and anothcr limiting case applies, namely SI ^> T. In such cases <|69|) the gravitational 
state decay within the lifetime: 

r = 2/r 

which is small compared to the passage time through the wave-guide. 

On FigEI we pl°t the results of numcrical calculations for the measured neutron flux within the time-dependent 
model for 2 vàlues of roughness amplitude 6=1 /im and b — 2 fim and | Imet| = 0.1 fim. Better resolution of the 
quantum "steps" appears with an increase of the roughness amplitude. 

In the above simple "two-state" model several potentially important effects are not taken into account, in particular 
the "non-resonant" transitions between different gravitational states. However this model enables understanding of fast 
irregularities ("steps") in the transmitted neutron flux as a function of absorber position H and naturally explains 
them in terms of gravitational states of neutrons. The model also establishes the dependence of the wave-guide 
absorbing properties on roughness amplitude. 

C. Resolution of gravitational states 

1. Constraints on resolution 

Based on the results of the previous sections we can analyze the conditions for the best resolution of gravitational 
states. The presented numerical calculations show that an increase in absorber efficiency (e.g. by increasing roughness 
amplitude) results in a shifting of the positions of "the quantum steps" in the neutron flux by the value A„ and 
enhancing their resolution S n . To perform a qualitative analysis we will accept that the "step-like" increase in the 
measured neutron flux, corresponding to the "appearance" of the new state, starts to be seen when the widths of this 
state are: 

T n (H)T pass =e 

We will also accept that such a "step-like" increase saturates when 

r n (H)rP ass = 1/e 

From expression (|70J) one obtains the following estimate for the shift A„ and the uncertainty S n of the n-th step in 
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the extreme limit \n(T pass /rf s ) > 1: 



A, 



2/3 



ln(rP as 7rf s VA^n) 



2/3 



3 14 



2/3 



ln(r pass /rf s VA^:) 



-1/3 



(71) 
(72) 



where 



Ao = /o(0 2/3 [ln(r^7r ï f s )] 2/3 



The uncertainty S n decreases as In -1 ^ 3 (T pass /T^ bs ) with an increase in T pass /t^ òs . The resolution of the n-th 
state is possible if the uncertainty in the step position S n is much less than the distance between neighboring steps 
H n+ i + A n+ i — H n — A„. For highly excited states we can use the WKB expression (JJJ) for the classical turning point: 



H n — lo 



3tt 



(271-1/2) 



2/3 



to find the universal limit on the number of states that can be resolved if T pass /t^ s > 1: 

[ln(T pafl Vr° 6 s )] ~ V3 < ^(2tt) 2 / 3 [(n + 3/4) 2 / 3 - (n - 1/4) 2 / 3 

This estimation shows that the resolution of states very slowly increases with an increase in passage time or in the 
efhcicncy of the absorber in the limit r pass /r^ bs » 1, namely n ~ ln(T pass /r° bs ). This law is the consequence of the 
linear dependence of the gravitational potential on z. Indeed due to the linearity of the gravitational potential, the 
level spacing decreases with n like n -1 / 3 , until the neighboring states' contribution to the flux starts to overlap. In 
particular for the value of T pass /t* = 100 the number of states that can be resolved is around 5. 

The resolution of quantum states could be improved, if the initial population of one or several of such states is 
artificially reduced. In this case the neighboring state would be exposed. We have studied the scenario in which the 
bottom mirror has a specially designed "step" 



íg state 

mm. 



2. Repopulation of states 

If two bottom mirrors are shifted relative to each other by a A of a few jj,m in height, there is an additional boundary 
at the step position x = Lo that will change the population of the eigenstates. We give here just a brief description 
(for the details see p0{). 



I X — Lq I X — Lq 



9 T 

ox 



x=L 



9 T 

ox 



Vz e [o, h] 



x=L 



Due to the presence of the shift A in the bottom mirrors' position, the gravitational states are repopulated. If this 
step is treated as a "sudden change" in the potential, the matching at the boundary x — Lo results in the following 
repopulation coefhcients: 



C jm = cxp (-r,-T /2)) / ipj(z)tp m (z + A)d 
Jo 



(73) 



Here t = Lo/V, <fj(z) is the gravitational state in the presence of the absorber, positioned at height H above the 
first mirror. Again, note the usage of the biorthogonality condition. 

The expression for the neutron flux at the detector position is now modihed as follows: 



(74) 



We neglect here the interference terms, assuming wide longitudinal velocities distribution. 



22 



n 




1 


0.162 


2 


0.765 


3 


0.037 


4 


0.019 


5 


0.009 


6 


0.005 


7 


0.002 



TABLE II: Normalized repopulation coeficients after transition from a single ground state across a mirror shift of m 9.15 /im 

To illustrate the effect of repopulation, consider a simplified system consisting of just two mirrors without any 
absorber. The orthonormal system of cigcnfunctions of the vertical motion in this case is just given by the Standard 
bound state Airy function. Now imagine the second mirror shifted downwards compared to the first one by an amount 
equal to the height of the first node of the 2 nd eigenstate wavefunction 

(A 2 — Ai) • lo ~ 1-56 • lo ~ 9.15 /ím 

It is clear that the 2 nd eigenstate above the 2 nd mirror exactly matches the ground state wave function above the l 3t 
mirror from its edge on, while the ground state wave function of the 2 nd mirror overlaps only with the exponentially 
decaying tail of the ground state of the l sí mirror. This implies immediately that the new ground state above the 2 nd 
mirror will be suppressed with respect to the 2 nd eigenstate above this mirror. The repopulation coefficients for the 
transition to the 2 nd mirror, normalized to the initial population of the ground state above the l st mirror, are given 
in TablcCH 

FigElcontains a plot of the population of the new ground state above the 2 nd as a function of the relative shift of 
the two mirrors. 

1 

0.8 
0.6 
0.4 
0.2 



2 4 6 8 10 12 14 
mirror shift [yim] 

FIG. 9: Repopulation coefficient of the ground state above the 2 nd mirror as a function of the relative shift of the two mirrors 
in fim 

In Fia llOl the neutron flux for mirror shift A = 8 \xm is compared with the neutron flux without any shift in the 
bottom mirror position. Due to the depopulation of the ground state, the changes in the flux slope corresponding to 
the gravitational states can more easily be seen. 

VI. CONCLUSIONS 

We have analyzed the problem of the passage of ultra-cold neutrons through an absorbing wave-guide in the 
presence of the Eart h's g ravitational field, both qualitatively and numerically. We have shown that the set of existing 
experimental results |l8l ll9l l20l | exhibits clear evidence for the quantum motion of neutrons in the gravitational field. 

We developed the formalism describing the loss mechanism of ultra-cold neutrons in the wave-guide with absorption. 
The essential role of the quantum reflection phenomenon for the loss of ultra-cold neutrons was established. The 
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FIG. 10: The relative neutron flux as a function of slit height for two vàlues of the bottom mirror shift. F* is the flux calculated 
at absorber position H = 35 jxm and diffuseness p = 1 /xm. 

concept of quantum reflectíem enables universal description of different kind of absorbers in terms of effective complex 
scattering length a. The efficiency of absorption of ultra-cold neutrons in the presence of the gravitational field of 
Earth is determined by the ratio of such a scattering length to the characteristic gravitational wave-length a/lo. 

We studied the particular case of absorbers with rough surface. It was established that in the latter case the main 
loss mechanism is duc to the non-specular reflection of neutrons from the rough edges of the absorber. Absorber 
cfhciency turns out to be proportional to the square of its roughness amplitude, if this amplitude is small compared 
to the characteristic gravitational wave-length Iq. 

We calculated the neutron flux through the wave-guide in the case of zero gravity (mirror and absorber arranged 
parallel to the gravitational field). For large slit hcights, the dependence of such a flux on the slit height H exhibits 
a power law. Its exponent depends on the absorber efficiency. These calculations are important for independent 
measurement of absorber/mirror properties. 

We argue the possibility of using the " inverse geometry" experiment for measuring the lifetime of neutrons bouncing 
on an absorbing surface. The neutron lifetime was found to be T abs = \/{2mg\ Ima|). It was determined by the 
gravitational force mg, acting on the neutron and imaginary part of the scattering length |Ima| of the absorbing 
surface Fermi potential. This experiment shows unambiguously the role of gravitation on the lifetime of ultra-cold 
neutrons. 

The theory developed in this paper allows to analyze the resolution of the gravitational spectrometer and to 
compare the efficiency of different kinds of absorbers/scatterers. We show that the spectrometer resolution is severely 
limited by a fundamcntal rcason: finitc penetrability of the gravitational barrier between the classically allowcd 
region and the scatterer height. The resolution can be improved by a significant increase in the time of storage of 
neutrons in quantum states, and/or by improvement of the efficiency of the absorber/scatterer. The efficiency of best 
absorbers/scatterers used in actual experiments was defined mainly by the shape of their rough surface so that the 
efficiency is approximately proportional to the square of the roughness amplitude (when the roughness amplitude 
is smaller than the characteristic scale of the gravitationally bound quantum states lo ). Further increase of the 
roughness does not improve the efficiency; however strict theoretical description of the case of a large amplitude 
roughness is not covered by the present analysis. Another way of increasing the resolution could be through the 
selective depopulation of certain gravitational states, for instance by applying a bottom mirror with a "step". 

The results obtained are rather general in character and can be applied to different physical problems, involving the 
transmission of quantum partides through absorbing wave-guides. The development of the theoretical considerations 
presented would include the incorporation of large roughness amplitudes comparable to or larger than the characteristic 
gravitational length lo <~ 6 /im; as well as the studies of long storage time case, when decay of neutron quasi-bound 
gravitational states differs from the exponential law. These are necessary if the highest resolution is to be achieved 
for the mcthod considered. 
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VIII. APPENDIX A 



Here we derive the equation for the energies of neutrons localized between an ideal mirror and absorber in the 
presence of a gravitational field. We assume that the absorber Fermi potential has a diffuse radius much smaller than 
the characteristic gravitational wave-length p <^l . In the region where the absorber potential can be fully neglected 
< z <C H — p the wave-function is the superposition of Airy functions: 

ip b (z) ~ Ai(z/Z - A„) - SBi(z/l Q - X n ) 

The zero boundary condition on the mirror gives: 

Ai(-À n ) 

S = TW \ X ( 75 ) 

Bi(-A„) 

The neutron wave-function inside absorber z > H — p is determined by the absorber Fermi potential, which is much 
stronger than the gravitational potential. In the range of distances H — Iq <ti z ^ H — p such a wave-function is 
weakly perturbed by gravitation and can be written as: 

, , x H-z 

Mz) ~ i + — — 



where à = a — H , with a being the complex scattering length on the absorber Fermi potential. One can see that à 
plays the role of the "scattering length of the diffuse tail" of the Fermi potential. 

Now we match the wave- functions ipb(z) and ip a (z) and their derivatives in the region H — lo<^z<tiH — p. For 
this we use the Tailor expansion of ipb(z) in the vicinity of H: 

Mz) ~ M(H/l - X n ) - SBi(H/lo - Xn) + {AÏ(H/l - X n ) -SBi'(H/l - X n )) [z - H) 

The matching condition gives: 

a = Ai(ff/Z - X n ) - à/l mM/h - A n ) , g] 

Bi(H/l -X n )-à/l Bi'(H/lo-X n ) 1 ' 



Putting together (|75(l and l|76(l we finally get the equation for the eigenvalues X n 

Ai(-A ?I ) = M{H/l Q - X n ) - à/l Ai' (H/lp - K) 
Bi(-A„) Bi(íf// - A„) - à/l BiWo - A„) 



(77) 
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The equation for the eigenvalues in an inverse geometry experiment can be obtained in a similar way. The wave- 
function outside the absorber ipb(z) now vanishes at the mirror position H, which gives for Si nv : 

_ M(H/l - A„) 

bmv ~ Bi(H/l -\ n ) U8j 
The wave-function ip a {z) of the neutron inside the absorber at the asymptotic distances z 3> p is: 

ip a (z) ~ 1 - z/a 

The matching of ^> a {z) and ipb(z) at distances Iq S> z 3> p together with l|78[l results in the following equation for A„: 

a [M{H/l - A n ) Bi'(-A n ) - Ai'(-A n ) Bi(H/l - A n )] = Ai(-A n ) Bi(ff/Z - A n ) - Bi(-A n ) Ai(H/l - A n ) (79) 

Note that the derivation of the above equations is based on the fact that p <C Iq, so that the wave-function tpa(z) 
is weakly perturbed by the gravitational field in the asymptotic region p <C z <C ^o- 

IX. APPENDIX B 

In this Appendix we derive the useful relation between the nonadiabatic coupling matrix element (Vjl^^) an d the 
energies of corresponding states i and j. 

We will study the one dimensional Schròdinger equation: 

01^)=^) (80) 

The eigenfunctions <fj{x) and ipi(x) obey the following boundary condition: 

<Pi(x = 0) = (81) 
l Pi {x = H) = (82) 

Here the varying parameter H is a boundary. Hereafter we assume that the Hamiltonian itself is independent of H, 
while eigenfunctions pi(x,H) and energy eigenvalues Ei{H) depend on H through the boundary condition (|82|l only. 
Applying d/dH to both sides of l|5ü|)we get: 

Integrating the left side of (|83[1 with (pj(x,H) and taking into account boundary conditions l|81l) and <|82|) we get: 

/ iúi d ^\ d(p j {x)d(p l difi - 



Note, that: 



Combining the above results we get for the matrix element of interest: 

{vj 1 1|) = * BH^B+ eH ò v (84) 

Now let us use the following relation 



dH 

From we get in case i = j: 



= 2(dipi/dH\(pi) = 



d(pi(x) difi(x) dEi 

~Tx — dïT [x=H --dH (85) 
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It is clear that the expression (|84|l can be expressed as: 



<VJI 



tjtj 



dH' 



Ej — Ei 



From (|85|l we finally get: 



(Vil— ) 



y/dEi/dHdEj/dH - dE t /dHS t 
Ej — Ei 



Applying the above result to the coupling matrix element in the time-dependent model l|59|) we get: 



a(H) = 



y/d\ n /dHdX*/dH 

Ari — A* 



(86) 



(87) 



Here X n is the eigenvalue of the low-lying gravitational state, while A* is the eigenvalue of the highly excited "box- 
like" state. This expression is much more convenient for practical applications than the integral in the definition of 
the coupling matrix element. In particular, it can be used to obtain the asymptotic expressions for the width (|7U|) of 
a given gravitational state n if H ^> H n . 

To obtain such an expression we first find the eigenvalue derivative d\ n /dH from the equation (|63l) : 



dX r , 



Ai' (H/l Q - A„) Bi(-A„) - BÏ(H/l - A„) Ai(-A„) 



8H l kí'{H/l - A„) Bi(-An) - BV(H/l Q - A„) Ai(-A„) + Ai(H/l - A n ) Bi'(-A n ) - Bi(H/l - X n ) Ai'(-A„) 

(88) 

Taking into account the expression l|63|l and asymptotic properties of the Airy function of large argument H/lo 3> A„ 
we get: 



dX n 
dH 



1 H-H n 



■ exp 



-4/3(JÏ/Z - A„) 3/2 



V H n 

For the energy E* of the highly excited "box-like" state, we can use expression (|41|l . from which we get: 



dE* 



E* 



dH H 

For the square of the coupling matrix element a 2 (H) in case of large H H n we get so far: 

E* so \H — H n 



arÇH) = 2 



H(E* - E n y lo V H n 



■ exp 



-4/3(íf/í - A„) 3 / 2 



Taking into account the expression for the width T* of the "box-like" state (|42|l and substituting the above results 
into the expression for the width of gravitational state (|68|l : 



r n = &Va 2 (ff)/(4r*) 



we finally come to the expression: 



r n — So 



H n 8lo \ Imal 



H - H n 



lo 



■ exp 



--{(H-H n )/lof 2 



